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Abstract 

A method of integrable discretization of the Liouville type nonlinear partial differential 
equations is suggested based on integrals. New examples of discrete Liouville type models 
are presented. 

1 Introduction 



H ' The problem of integrable discretization of the integrable PDE is very complicated and not enough 
studied. The same is true for evaluating the continuum limit for discrete models [T]. In the 
present paper we undertake an attempt to clarify the connection between Liouville type partial 
differential equations and their discrete analogues. One unexpected observation is that there are 
pairs of equations, one continuous and the other one semi-discrete, having a common integral. 
Inspired by these examples, we introduced a method of discretization of PDE having a nontrivial 
integral. Similar ideas are used in [2] where a method of construction of difference scheme for 
ordinary differential equations preserving the classical Lie group is suggested. Let us begin with 
the necessary definitions. 



^e-mail: habibullinismagil@ginail.com 
^ e-mail : natalya@fen. bilkent .edu.tr 
■^e-mail: alfiya85.85@mail.ru 



We consider discrete equations of the form 



v{n + 1, m + 1) = f{v{n, m),v{n + 1, m),v{n, m + 1)) (1) 
and semi-discrete chains 

t{n + 1, x) = /(x, x),t{n + 1, x),tx{n, x)) . (2) 

Equations ([T]) and ([2]) are discrete and semi-discrete analogues of hyperbohc equations 

Uxy = f{x,y,u,u^,Uy) . (3) 

Functions v = v{n,m), t = t{n,x) and u = u{x,y) depend on discrete variables n and m and 
continuous variables x and y. Through the paper we use the following notations: 

Vij = v{n + i,m + j); Vi = Vi^; Vj = Voj] ti = t{n + i,x). 

For equation ([3]) function W{x,y,u,Uy,Uyy, . . . ,d^u/dy'') is called an x-integral of order k if 
D^W = and Wgh^/dyk ^ 0, and function W{x, y, u, u^, u^x, • • • , d"^u/dx"^) is called a y-integral 
of order m if DyW = and Wdmu/dx"^ 7^ 0. Here, and Dy denote the total derivatives with 
respect to x and y. Equation ((31) is called Darboux integrable if it possesses nontrivial x- and y- 
integrals. 

For equation ([2]) function F{x,n,tm,tm+i,tm+2, ■ ■ ■ ,tm') is called an x-integral of order m' — 
m + 1 if Dj.F = and Ft^ ^ 0, F^/^ ^ and function I{x,n,t,tx,txx, • • • , ^) is called an n- 
integral of order k, if DI = I and I^kt 7^ 0. Here, D is the forward shift operator in n, i.e. 
Dh{n, x) = h{n + l,x). Equation ([2]) is called Darboux integrable if it possesses nontrivial x- and 
n- integrals. 

For equation ([T]) function /(ra, m, Vk-, Vk+i, . . . ,Vk') is called an n-integral of order /c' — /c + 1 if 
DI = I and ly,^ 7^ 0, ly^, 7^ 0, and function J(n, m, f^, f^+i, fr+2, • • • , "^r') is called an m-integral of 
order r' — r + 1, if /)/ = / and 1^^ 7^ 0, 7^, , 7^ 0. Here, D and Z) are the forward shift operators in 
n and m respectively. Equation ([T]) is called Darboux integrable, if it possesses nontrivial n- and 
m- integrals (see also [3]). 

Continuous equations ([3]) are very-well studied. In particular, the question of describing all 
Darboux integrable equations is completely solved ( see |1] - [7]). All equations possessing 
X- and y-integrals of order 2 are described by the following theorem. 



Theorem 1.1 (see ^Tj) Any equation for which there exist second order x- and y-integrals, 
under the change of variables x — > X{x), y — > Y{y), u — i- U{x,y,u), can be reduced to one of the 
kind: 

(1) U^y = 1? = U^^ - Q.bul, W = Uyy- CSM^/ 

(2) u.,y = eyuy, W = u,-e^,W = ^- Uy; 

(3) u^y = e-^ul-A, W = Uxx- 0.5^2 - O.Se^", W = ^^^q^/ 

U) u^y = UxUy ^-], W = ^-^ + ^,W=^-^ + -^; 

\-t/ ^ y yu—x u—y J ' Ux u—x u—x' Uy u—y u—y' 

(5) u^y = ij{u)f3{u^)(3{uy), {Impy = (3(3' = -u^, (3(3' = ~Uy, 

(6) Uxy = PJ^^^Ail^ pp' + c(3 = -Ux, W + cP = -Uy, 

13 u' p u' 

(V Uxy = -2^, W=^ + 2^,W=^ + 2^^; 

(8) Uxy = ^x^y),lxmuyy P' = P' + /3^ P' = P' + ^^ 

W = U^^(3{U^) - J^^^^J^, W = UyyP{Uy) - (^^^J^ ■ 

On the contrary, the problem of describing all equations (II]) or ([2]) possessing both integrals (so- 
called Darboux integrable equations) is very far from being solved (the problem of classification is 
solved only for a very special kind of semi-discrete equations |8]), it would be beneficial for further 
classification to obtain new Darboux-integrable equations ([T]) and semi-discrete chains (E]). It was 
observed that many chains (|2]) and their continuum limit equations (|3]) possess the same n- and 
y-integrals: 



semi — discrete chain 


n — integral I 


continuous 
analogue 


y — integral W 


hx = tx + 0.5t? - 0.5^2 


t^ - O.St^ 


Uxy UUy 


Ux — O.Su^ 


t^^ = t^ + C'e0-5(*+*i), C = Const 


txx 


Uxy ^ 


u^x - 0.5ul 


hx = tx + Ve2* + Ce*+*i + e2ti 


t^^ - O.btl - O.Se^* 


Uxy = e"yi + 


u^x - O.bul - 0.5e2« 



The main aim of the present paper is the discretization of equations ([3]) preserving the structure 
of y-integrals of order 2: we take y-integral for each of eight classes of Theorem 11.11 and find the 
semi-discrete chain possessing the given n- integral (y- integral) . The next Theorem presents a 
list of semi-discrete models of Darboux integrable equations (|3]) from Theorem 11.11 with integrals 
of order 2. 



Theorem 1.2 Below is the list of equations ^ possessing the given n-integral I : 



given n — integral 


the corresponding chain 




I = t,.,. - 0.5^2 


t^^ = t^ + ^e°-5(*+*i), C = Const 


(1*) 


/ = - e* 


t\x ^x ^ ~h 6 


(2*a) 


J ^xx -f- 

I'X 


ti^ = K(t, ti%, where KtR-^ + Kt, = K - 1 


(2*6) 


I = t^^- 0.5tl - 0.5e2* 


tlx = t^ + Ve2* + Re^+ti + e2*i, R = Const 


(3*a) 




ti^ = (1 + Re'+''% + Vi?2e2(t+ti) + 2i?e*+*yt2 _ 4 


(3*6) 


— ixx ^ix 1 1 


- (t+L)(t-x) ^x, L - Const 


(4*) 


{Inipy = ij\f3f3' = -t. 


P{tx) = itx and tix = K{t,ti)tx, where 
Kt + KKt, + K^i^iti) - K^{t) = 


(5*) 


I3{tx) t ' 

+ c/3 = -t. 


(3{tx) = Rtx, and ti^ = K(t,ti)tx, where 


(6*) 


J ixx 1 '^y/ix. 


tlx = [y/t^ + ^^7^) , R = Const, C = Const 


(7*) 


I — l^{tx)txx {x+R)l3{ix)' 


/3{tx) = -1 and ti, = tx + 


(8*) 



It is remarkable that each equation in Theorem 11.21 also admits a nontrivial x-integral. It means 
that discretization preserving the structure of ^/-integrals sends Darboux integrable equations ([3]) 
into Darboux integrable chains ([2]). 

Note that equation (1*) was found in [9]. Equation (3*a) for i? = 2 was found in [3], equations 
(2*a) and (3*a) are found in |H]. To our knowledge, the other equations from Theorem 11.21 are 
new. 

The next theorem lists x-integrals for chains from Theorem II. 2[ 

Theorem 1.3 (I) The equations (2*6), (5*) and (6*) from Theorem having the form ti^ = 
K{t,ti)tx admit x-integral F{t,ti), where function F is a solution of Ft + K{t,ti)Ft^ = with a 
given function K{t,ti). 

(II) X -integrals of equations (8*), (1*), (3*a), (3*6), (4*), (7*) and{2*a) are F = {ti-t+C)/{x+y), 
p ^ g(ti-t)/2 _^ g(ti-t2)/2^ p ^ arcsinh{ae^^^^-^ + 6) + arcsinh{ae^^'^ + 6) with a = 2(4 - i?2)-V2^ 
6 = i?(4 - i?2)-i/2^ p = Vi?e2*i + 2e*i-* + Vi?e2*i + 2e*i-*2, F = (ti - t){t2 + L){t2 - t)-\ti + 
L)-\ F = {2ti - t - t2)/(2C2) - l/{x + R) and F = (e* - e'^){e'' - e*3)(e* - e*«)-i(e*i - e'^y^ 
correspondingly. 



One can also apply the discretization method preserving the structure of integrals for semi- 
discrete chains ([2]): take x-integral for a semi-discrete chain and find discrete equation ([1]) with 
the given m-integral (x- integral). 

In spite of the absence of the complete classification for Darboux-integrable semi-discrete chains 
([2]) there is a large variety of such chains in literature (see, for instance, and [IQ]). The 

procedure of obtaining fully discrete equations for a given integral is a difficult task and requires 
further investigation. As a rule it is reduced to a very complicated functional equation. We 
illustrate the application of the discretization method on chains (1*), (4*) and (7*) from Theorem 
11.21 The discrete analogues of the chains are presented in the next Remark. 



Remark 1.4 Below is the list of equations (Cp possessing the given m-integral I : 



given m — integral 


the corresponding equation 




I = e(^i"^)/2 -I- e(fi-'"2)/2 


I-.V-I.T+V _ 1 


{!**) 


I = (vi- v){v2 + L){v2 - v)-^{vi + L)-^ 


L{vi+vi-v)+vivi 
^1.1 - L+v 




I = 2Vi — V — V2 


Vi^i = f 1 + h{vi — v), z = h{2z — h{z)) 


(7**) 



The equations {!**), (4**) and (7**) have respectively the following n-integrals I = e^'"'^ ^^1"^ + 
g(0i-i;2)/2^ / = - v){v2 + L){v2 - ty^{vx + L)'^ and I = vx-v - h'~^{vx - v) with h'^ being 
the inverse function of function h that satisfies the functional equation z = h{2z — h{z)). 



Equation (1**) from Remark 11.41 appeared in [IT], equations (4**) and (7**) seem to be new, 
unfortunately we failed to answer the question whether equation z = h{2z — h{z)) has any solution 
different from linear one h{z) = z + C. 

The article is organized as follows. Theorem 11.21 is proved in Section 2. The proof of Theorem 
11.31 is omitted. Chains (1*), (2*a) and (3*a) are of the form tix = tx + d{t,ti), and their x- 
integrals can be seen in [B]. One can find x-integrals for chains (3*&), (4*), (7*) and (8*) by direct 
calculations. In Section 3 the discretization of chains (1*), (4*) and (7*) from Remark 11.41 are 
presented and for each obtained discrete equation the second integral is found. In Section 4 the 
Conclusion is drawn. 



2 Proof of Theorem 



1.2 



Case (1*): Consider all chains ([2]) with n-integral of the form I = t^x — V^x^- Equality DI = I 



implies 



fx + fttx + ftif + ft^txx — -^f"^ — txx — ■ (4) 



By comparing the coefficients before t^x in © we have ft^ = 1. Therefore, 

f{x,t,ti,tx) =tx + d{x,t,ti). (5) 

We substitute (E]) into (jl]) and get dx + dttx + dt^tx + dt-^d — ^tx'^ — dtx — ^cP = —\tx'^, or equivalently, 
dt + dt^ — d = and dx + dt^d — ^d"^ = 0. We solve the last two equations simultaneously and find 
that d = e^^K(x,ti — t), where K = Ce^^*^*^^*-* and C is an arbitrary constant. Therefore, chain 
([2]) with n-integral I = txx — \tx^ becomes tix = tx + Ce(*i+*)/^. 

Case (2*a): Consider all chains ([2]) with n-integral I = tx — e*. Equality DI = I implies / — e*^ = 
tx — e^, which gives the equation tix = f = tx — + 6*^. 

Case (2*6): Consider all chains (|2]) with n-integral I = — tx- Equality DI = I implies 

fx + fttx + /ii/ + ft^ txx f txx . /^\ 

7 J = -. tx- (6) 

/ '^x 

By comparing the coefficients before t^x in dH]) we have ft^/f = 1/tx, that is / = K{x,t,ti)tx- 
Substitute / = K{x, t, ti)tx into (E]) and have ^ + ^tx + Kt^tx — Ktx = —tx, or equivalently (by 
comparing the coefficients before tx and tx^), we get + Kt-^ = K — 1 and Kx = 0. Therefore, 
equations tix = Kit, ti)tx, where K satisfies + K^^ = K — 1 are the only chains ([2]) that admit 
n-integral / of the form / = ^ — tx- 

Case (3*a): Consider all chains (j2]) with n-integral I = txx — — 1^^*- Equality DI = I implies 

fx + fttx + ftif + ft^txx — -/^ — -e'^*^ = txx — -^tx"^ — -e^*. (7) 

By comparing the coefficients before txx in dB we have ft^ = 1, that is f{x, t, ti, tx) = tx+d(x, t, ti). 
Substitute f{x, t, ti, tx) = tx + d{x, t, ti) into ([7]) and have 

dx + dttx + dt.itx + d)-]^{tx + d)^ - ^e^*^ = -^t,.' - ^e'*. (8) 

Compare the coefficients before tx and tx^ in ^ and get 

dt + dt,-d = Q, dx + dt,d-^d^ -^e^'' = -^e^'. (9) 



The first equation in (Q has a solution d = e^^K{x,ti — t). Substitution of this expression into 
the second equation of (Q gives e'^^K^ + Kt^_tK + ~ \ + |e"^(*i"*) = 0. Since K depends 
on t/ = ti — t and x, then = and the last equation becomes 2K'K + = 1 — e~^^, and 



hence, d = e^^K = a/ e^*! + e^* + i?e*+*i , where is and arbitrary constant. Therefore, chain ([2]) 



with n-integral I = t^x — \tx^ — ^e^* becomes tix = + V e^*! + e^* + i?e*+*i , i? = const. 
Case (3*6): Consider all chains (|2]) with n-integral I = h^^^^^. Equality DI = I implies 



By comparing the coefficients before txx in ( |TOl) we get ^y^— = — ^===, that is arccosh^ 
arccosh^ + /^(x, t, ti). Thus, 



/(x,t,ti,t,) = At^ + ^Vtx -4, (11) 
where A(x,t,ti) = coshK, B{x,t,ti) = sinhK, — B'^ = 1. Note that / = 2cosh{{arccoshif) + 



/r), i.e. VP -4^ = 2sinh{{arccoshif) + K) = 2{J^-f- - IcoshK + ifsinhK), or V/^ - 4 



Et^ + A _ 4 Substitute (III]) into ([TO]) and have 



+ Bx^tx^ - 4 + + t.^tV^x' - 4 + {txAt, + Bt.^tx^ - 4) (At, + B^tx^ - 4) 

-{Atx + - 4)2 + 4 = -(5t,. + Av/t.2-4)v't.'-4, 

that can be written shortly as 

{tx^ - 4)(ai + a2txf = {as + a^tx + a^tx^f , (12) 

where ai = Bx, a2 = Bt + At^B + Bt^A - 2AB + B,as = -ABt^B + 45^ + 4 - 4^4, 04 = Ax, = 
At + At^A + Bt^B - A^ - B"^ + A. We compare the coefficients before tx^, t,^, tx, in 
and have = ^5^, 2aia2 = 20405, — 4q;22 = a^^ + 20305, —80102 = 203O4, —401^ = 03^, 
that implies Oi = 02 = 03 = 04 = 05 = 0, which is possible only if A = 1 + i^e*"*"*^ and 



B = -\/_R2g2(t+ti) _|_ 2i?e(*+*i), where R = const. Therefore, by (ITT]) , the chain 02]) with n-integral 



/ = +^ becomes ti, = (1 + /2e*+*i)t, + Vi?2g2(t+ti) + 2i?e(*+*i)^t,2 _ 4_ 

Case (4*): Consider chains (E]) with n-integral / = 7^ - ^ + Equality DI = I implies 

fx + /ttj;.. + ftlf + ft^txx _ 2f ^ 1 _ txx^ ^ 2tx ^ 1 /-j^gN 

/ ti— a; ti — X tx t — X t — X 

We compare the coefficients before txx and have ft^j f = l/tx, that is / = txK{x,t,ti). Substitute 
f = txK into ( IT3]) and have 

+ i^tt,' + i^fii^t,' 2irt, ^ 1 _ 2t, ^ 1 

i^'tr t-t — X t-i — X t — X t — X 



By comparing the coefficients before tx and t^^ in (fT4|) we get 

Kt _ 2K 2 Kx _ 1 1 

Y ~ ~ 'k ~~t^^T^' ^^^^ 

We solve two equations of f|T5|) simultaneously and have K = ^^^f^, where L is an arbitrary 

- + ^ b 



constant. Therefore, any chain ([2]) with n-integral / = ^ — ^i^- + -1- becomes ti^ = ^jTrH~^^ 



Case (5*) : Consider all chains ([2]) with n-integral I = ^ where /3 = /3{tx),ip = il){t),f3f3' 



-tx. We have, 2/3/3' = -2t^, i.e. /32 = -t^^ + M^, or /3 = - where M is an arbitrary 

constant. Equality DI = I implies 

^ ^^^^ ^(^^^^^ ^ - v^(ti)/3(/) = ^ - mm- (16) 

We compare the coefficients before t^x and have ft^/f3{f) = l/f3(tx) which implies that 
either (5*a): M = 0, /3()f:i.) = itx and ti^, = K{x,t,ti)tx, 
or (5*6): M 7^ and then arcsin-fr = arcsin¥j + t,ti), that is. 



/ = txA{x, t, ti) + - tx''B{x, t, ti), ^2 + fi2 = 1. (17) 

In case (5*a) we substitute tix = f = K{x,t,ti)tx into f lT6|) . use that /3(ta,) = it^., and obtain 

= 0, ^ + Kt,+^{h)K = 4j{t). (18) 

Therefore, the chains ([2]) with n-integral / = — iip{t)tx are equations ti^; = K{t,ti)tx, where 

function K satisfies ( !T8|) . 

Let us consider case (5*6). Note that 



-f = M^- AHx^ - 2ABtxJM^ - tx^ - B^M^ + BX^ = {Btx - AJm^ - t^2)2 



and /3(/) = ±{Btx - AJM^ - t^.^), /3(t^) = a/M2 - t^.^. Substitute ([IT]) into (116]) and get 



Axtx + BxJM^ - tx^ + Attx" + BttxJM^ - tx" + {At.tx + 5*, JM^ - t,2)(^^^ + qJ^^ _ ^ 



fx J 



±(5t, - Aa/M^ - tx") 



= ±{Btx - A^M^ - tx^)ij{h) - - txXit), 

or the same, 

(M^ - t,2)(ai + a2t^,)2 = (^3 + a^tx + a^tx'f, (19) 

where ai = Bx, = Ax, aa = Bt+At,B+ABt^+2ABilj{ti)+Bip{t), = BBt^M'^-A^M'^%lj{ti)- 
Atlj{t)]VP, a5 = At + At^A - Bt^B - B'^tpiti) + A'^ipiti) + Aipit). We compare the coefficients 



before tj;, k = 0,1,2,3,4, in (fT9|) and find that ai = ^2 = «3 = «4 = «5 = 0, which is possible 
only if = i? is a constant function, that contradicts to the equation {Inip)" = Therefore, 
case (5*6) is not realized. 

Case (6*): Consider chains with n-integrals / = - where /3 = (3{Q and (3(3' + c(3 = 



—tx- The equality DI = I implies 

ft:c 



and 

f. + ftt. + ftj _PU)_^ Pit.) 
(3{f) h t ' 

Differentiation of ( 120|) with respect to x, t, ti gives 



(20) 



(21) 



First we differentiate (!2T|) with respect to t^, use (122|) . and get 

1 1 ;,.._M|)±Z1^£^^. p3) 



Next we differentiate ( l23ll with respect to t^;, use (l22l) . and arrive to the equality 



t. f \ f _ + f% m ^ Pit.) ^ ^ ti 

jti — , n/. \ : — ' — : — ^ — + 



[Pit.) Pif)}-" ti/3(y ti t t tpit^y 

There are two possibilities: 



either (6* a), when 



or (6*6), when 



mpit.) hicPif) + f)t. Pif) Pit.) ti \ 
t.Pif)-fPit.)[ t, t ^ t ^ ' 

Let us first consider case (6*a). It follows from (!24|) and ( 120|) that ft^/f = l/t., that is / = 
K{x,t,h%. We substitute / = K{x,t,ti)t, into (EH), use (3{f)/ti = iPiQf)/itJi) = 
and obtain 

that is, = 0, (3{tx) = \J R^tl + Ctx-, R = Const, B = Const, and 



Substitution of (3{tx) = \J B?t1 + Ctr^ into shows that (3{tx) = Rtx- Therefore, in case (6*a), 
the n-integral is / = ^ — ^ and the corresponding chain ([2]) is of the form ti^ = K(t,ti)tx, 
where K satisfies f l26|l . 

Let us now study case (6*6). It follows from ( l25l) and ( l23l) that 



First we differentiate (!25|) with respect to t and find /^^j, use the expression for ft from (j2j 
and /?'(/) = -(/ + c(5{f))/ I3{f) to express /tji in terms of /?(/), /^(t^:), /, t, ti, t^- Then we 
differentiate (1271) with respect to ti and find /^t^, use the expression for ft^ from (l25l) and /?'(/) = 
-(/ + cl3{f))/l3{f) to express /ti^ in terms of /3(t^), /, t, ti, t^. 
Direct calculations show that 

. _ 2(3if)cW{f) + + f')i-tf + hicPjtx) + tx)) 

Equality /t^, = ft,t yields (i) + c//3(/) + = 0, i.e. /?(/) = Af , f3{tx) = Atx, where 



A = or (ii) / = ht~\c/3{tx) + tx). 

Let us consider case (i). It follows from ( !20l) that / = K{x,t,ti)tx. The same considerations 
as in part (6*a) show that the chain ([21) in this case is tix = K{t,ti)tx, where function K{t,ti) 
satisfies fl26|) . 

Let us consider case (ii). It follows from ([20]) that /3(/) = ht-^{{l - c^)l3{tx) - ctx). We 
substitute this expression for /?(/) into (ED) and get c'^{2-c^)l3^{tx) + 2c{l-c^)txl3{tx)-cHl = 0, 
that implies that (I) c = 0, (II) = 2, (III) /3(t^.) = ^tx, or (IV) /3(tx0 = -^^x- Cases (II) and 
(IV) are not realized, each of them is incompatible with /3/3' + cP = —tx. Case (III) is realized only 
for c = 2 (with I3{tx) = —tx) and c = — 2 (with f3(tx) = tx). Therefore, using / = tit~^{cf3(tx) +tx) 
and the fact that c = (with (3{tx) = ^itx) or c = ±2{l3{tx) = — ± t^) we arrive to a chain ([2]) 
of the form tix = iy^x- Note that chains tix = ±tit^^tx with (3{tx) = ±t^, or P{tx) = ^itx is of 
the form tix = K{t,ti)tx, where K satisfies fl26l) with R"^ = 1 (for tix = —tit^^tx) or = —1 (for 

tlx = tlt~^tx)- 

Case (7*): Consider chains (12]) with n-integral J = ^ + 2^, y = Const. Equality DI = I 



implies 



fx + fttx + /ti/ + ft^txx ^ 2 ^xx ^ 2 (28) 



By comparing the coefficients before txx we have ft^j \pf = 1/y/t^, or 



f = {Vt:+K{x,t,ti))^. (29) 



Substitute ([29]) into ([28]) and get + KtU + KuU + 2Ku^K + K^K^ + ^ = 0. We compare 

x~\-y 

the coefficients before -y/t^, t^, tj^ and have 2Kt^K = 0, i.e. K = L{x,t); Kt + Kt^ = 0, i.e. 
K = L{x); and + Kt,K^ + ^f^ = 0> i-e- ^ = = const. Therefore, chain with 

n-integral / = + 2^^ becomes tix = {y/t^ + ;^)^5 where C and ?/ are arbitrary constants. 
Case (8*): Consider chains ([2]) with n-integral I = (5{tx)txx — {^x+y)i3{t ) ' where y is an arbitrary 
constant and P'itx) = f3^(tx) + The equahty Z^J = / gives 

+ /ttx. + ftj + ftj..} - , , = - ^ 



that imphes 

m)h. = Pit.) (30) 

and 

Pi,mf^,f,^^u}^^-^^-^-^^. (31) 

Differentiate (l30i) with respect to x, t, ti and get 

/x*. = -(/?(/) + /it. = + /m. = -(/?(/) + (32) 

Now differentiate (13T1) with respect to t^, we have 

a; + y (x + ?/)/^(/) 

Differentiate fl33]) with respect to ta, and get /t^ = — i^x+y)i3{f) ■ '^^^ ^^^^ equation together with 
([33]), ([30]) and ([31]) gives 

f = ^ f = 1 f = ^ 

ix + y)m' ix + y)Pify /?(/) ^ ' 

and 

= ^ { W) ~ Kf)m) 'W)^ W)} ■ ^^^^ 

Since, by and ([35]), ft^x-fxti = ■^^jj^{f3{f)+l), then /?(/) = -1, and, therefore, by we 
have = (x+y)-'^, ft = -(x+?/)"\ ft^ = 1. Hence, f{x,t,ti,tx) = tx + ^+C{x). We substitute 
this expression for / into (135]) and obtain C{x) = C{x + y)^^, where C is an arbitrary constant. 
Therefore, with the n-integral I = + 2^ the chain becomes ti^ = tx + ^ + C{x + y)~^, 
where y is arbitrary constant. 



3 Proof of Remark 



Case 1** : Consider all equations ([T]) with m-integral / = e^^~^ + e^^~^^. Denote by e~^^ = Wj, j = 
0, 1,2, and e~'"^ = wi. In new variables / = ^^^^ is an m-integral of equation wi^i = g{w,WiWi). 
DI = I implies 

W2 + W gi + wi 

= • (36) 

wi g 

We differentiate both sides of fl5B]) with respect to W2 and apply the shift operator D^^, we have 



Therefore, 



9wi 

wi g \wiJ \ 9 w 



9 = hc(w,wi), gi = Vc[wi,g). (37) 

W Wi 



We substitute f l37|) into f p6|) and get 



w 



9 — = c{wi,g) + wi. (3J 

Wi 



Substitution of f p7|) into fl5S]l implies that 0(1(7,^1)11; = c{wi,g)wi, or the same, c{w,Wi)w = 
D{c{w,Wi)w). Suppose that equation i = g{w,WiWi) does not admit an m-integral of the 
first order, then c{w,Wi)w = D{c{w,Wi)w) = C = const. Thus, c{w,Wi) = C/w. Finally, 
g{w,wi,wi) = ^^^^ + Cw~^. Therefore, the equations ([1]) with m-integral / = e^^""" + e"'^~"^ 
becomes e^^-'^'^'" = (C + e~^'"'^^'"'^^)~^ , where C is an arbitrary constant. Note that this equation 
is symmetric with respect to variables vi and vi. Therefore, n- integral for the equation can be 
obtained by simply changing in m-integral variables vj into variables Vj, j = 1,2. 
Case 4** : Consider equations ([I]) with m-integral I = . Equation vi^i = f{v,Vi,Vi) can 

be rewritten as f = r{v, t>_i, vi). Equality DI = I implies 

f-v,)(f, + L) ^ (v,-v)(v, + L) 

{f,-v,){f + L) {v2-v){v, + Ly ^ ^ 

Take the logarithmic derivative of ( 139|) with respect to V2 and then apply the shift operator -D^^, 
we get 



(40) 



h + L h-v, V2 + L V2-V (f + L){f-r) {v, + L){vi - v^i)' 

We conclude from the second equation of (HUj) that 

f + L vi + L _ 

= K{v,vi}. 41 

/ - r vi-v^i 



Take the logarithmic derivative of ( HTj) with respect to f_i and get f — r = r„_j(fi — f_i). 
Differentiation of the last equality with respect to vi yields f^^ = r^_-^. We differentiate fHOj) with 
respect to v_i and use the fact that f^-^ = r„_-^, we obtain /^^ = ± ^^Jj ^ . 

First assume that fy-^ = — ^ . We have, f — r = D{v,V-i,vi){vi — v^i)^^. It follows from 
r^, -^ = — ^ that f—r = C{v, t>i, Vi){vi—v_i)~^ , and, therefore, f—r = C{v, Vi){vi—v^i)^^. We 
substitute this expression for f—r into f l4T|) and see that /+L = C(f , i;i)fC(t>, ■Ui)(t>i+L)(t>i— t>_i)~^ 
which is impossible since / does not depend on 

Now consider the case when fy^ = ^JjJ ^ . We have, f — r = {vi — V-i)D{v,V-i,vi). Also, 
fv-i = yfSj ^ implies that f — r = (f i — v_i)C{v,Vi,vi). One can see that D{v,v_i,Vi) = 
C{v,Vi,Vi) =: C{v,Vi). Therefore, f — r = C{v,Vi){vi — f-i). It follows from f HT]) that 

f = A{v,Vi)vi + Aiv,v,)L-L, (42) 

where A = CK. Note that A = A{v,Vi) and Ai = A{vi, f{v,Vi,Vi)). Substitute into ([31]), 
get 

{Avi + AL-L- vi){AiV2 + AiL){v2 - v){vi + L) 
= {A1V2 + All - L- vi)iAvi + AL){vi - v){v2 + L), 
and compare the coefficients before f|, we have 

Ai{Avi + AL-L- vi){vi + L) = Ai{Avi + A){vi - v). (43) 

It follows from (H3!) that Ai = or, by comparing the coefficients before Vi, one gets A = 
Therefore, by ( |42|) . we have the equation fi^i = / = ^('"'^+'"^'^^+^^'"^ _ Note that the equation is 
symmetric with respect to variables Vi and Vi. This observation allows one to write down an 
n- integral J by a given m- integral I by changing in I variables Vj into variables Vj, j = 1, 2. 
Case 7** : Consider all equations ([1]) with m-integral F = 2v — Vi — f_i = D^^I, where / = 
2vi — V — V2. Equation t>i i = /(t>, t>i, -Ui) can be rewritten as = r(t>, t>_i, -yi). Equality 

DF = F implies 

2tyi - f - r = 2v - vi - v^i. (44) 

We apply ^ and to (jHj) and find that f^-^ = 1 and r„_^ = 1. Therefore, f = vi + h{v,Vi) 
and r = t>_i + q{v, Vi). Substitute these expressions for / and r into (jS]) and get 



q = 2v\ — 2v — h. 



(45) 



Equation vi^i 



f = Vi + h{v , Vi) can be rewritten as 



vi = V + h{v-i,v-i^i) = V + h{v-i,v-i + q{v, vi)). 



(46) 



First differentiate (H6ll witli respect to f_i and tlien apply tlie sliift operator D^^, we get D~^hy + 
D^^hy-^ = 0, tliat is h = h{vi —v). Equations (l44l) - (H6l) give -Ui — f = /i(2'i;i — 2v — h), or by taking 
e = tJi — f one gets e = /i(2e — h{e)). Tlierefore, tlie equation witli m- integral / = 2vi — v — V2 
becomes fi.i = vi + h{vi —v), where h solves a functional equation e = h{2e — h{e)). This equation 
^^1,1 = Vi + h{vi — v) admits also an n-integral. Since the equation is of the form Dz = h{z) with 
z = Vi — Vi then we have D{z — h^^{z)) = z — h^^{z). Actually, D{z — h^^i^z)) = D{z) — z = 
h{z) — z = z — h~^{z) = z — h^^{z). Here we use the identity h{z) — z = z — h^^{z) which is 
equivalent to the functional equation z = h{2z — h{z)). 

4 Conclusions 

The problem of discretization of Liouville type equations is discussed. Besides purely theoretical 
interest as a bridge between two parallel realizations of the integr ability theory, this subject has an 
important practical significance. There are two-dimensional Toda field equations corresponding to 
each semisimple or of Kac- Moody type Lie algebra (see [12], [13] )• The question is open whether 
there exist integrable discrete versions of these. Different particular cases are studied in [H], |15j . 
[T6] . In the article a step is done towards the solution of the problem. An effective method of 
discretization is suggested based on integrals. It is known that the Backlund transform is a kind of 
discretization (see [3], [H]). We would like to stress that our method of discretization essentially 
differs from that one. Even though for some exceptional cases the semi-discrete equation obtained 
realizes the Backlund transformation of the original equation for the other examples it is not the 
case. 
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